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Abstract 

Introduction:  In  this  study  we  explored  direet  and  indireet  measures  of  autonomie  nervous  system  fimetion,  as  well  as  ehanges  in 
eardiovaseular  eomplexity,  during  hemorrhagie  shoek  (HS). 

Methods:  HS  was  indueed  in  anesthetized  sheep  (n  =  8)  by  removing  40  ml/kg  of  blood  in  four  10  ml/kg  steps  over  40  min.  Resuseitation  was 
performed  with  laetated  Ringer’s  and  re-infusion  of  shed  blood.  Renal  sympathetie  nerve  aetivity  (RSNA)  was  measured  by  mieroneurography. 
Speetral  analysis  of  heart  rate  variability  (HRV)  employed  fast-Fourier  transformation  of  the  R-to-R  interval  (RRI)  of  the  EKG.  This  generated  the 
normalized  high-frequeney  (HFnu)  and  low-frequeney  (LFnu)  powers  of  the  RRI,  and  their  ratio  (LFnu/HFnu,  a  proposed  index  of  sympatho-vagal 
balanee).  Additionally,  non-linear  methods  were  applied:  RRI  eomplexity  was  measured  by  approximate  (ApEn)  and  sample  (SampEn)  entropy 
methods;  RRI  fraetal  dimension  was  measured  by  eurve  lengths  (FDCL).  Plasma  eateeholamines  were  determined  by  HPLC. 

Results:  The  model  eaused  profound  HS;  2/8  animals  survived  till  the  end  of  resuseitation.  RSNA  inereased  in  7/8  sheep  and,  as  HS 
progressed,  multiple  burst  eomplexes  were  identified  followed  by  sympathetie  withdrawal.  Coneomitant  deereases  in  HFnu  and  inereases  in 
LFnu/HFnu  oeeurred  after  20  ml/kg  blood  loss.  ApEn  and  FDCL  deereased  after  withdrawal  of  40  ml/kg  of  blood.  Cateeholamine 
eoneentrations  inereased  throughout  HS.  LFnu/HFnu  and  RSNA  were  not  linearly  eorrelated. 

Conclusions:  HS  led  to  an  inerease  in  RSNA  with  subsequent  withdrawal.  LFnu/HFnu  inereased  during  HS  in  assoeiation  with  vagal 
withdrawal  and  loss  of  RRI  eomplexity.  RRI  eomplexity  may  in  part  refleet  vagal  modulation  of  the  heart  rate.  Changes  in  direetly  measured 
tonie  sympathetie  traffie  do  not  eorrelate  with  non-invasive  measures  of  autonomie  modulation  of  the  heart. 
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1.  Introduction 

The  autonomic  response  to  hemorrhagic  shock  (HS)  is 
complex,  under-investigated  and  influenced  by  the  dynamics 
and  magnitude  of  volume  depletion  (Secher  et  al,  1 992;  Evans 
et  ah,  2001).  Animal  studies  involving  direct  renal  neural 
recordings  during  HS  showed  increased  renal  sympathetic 
nerve  activity  (RSNA)  as  the  initial  response  to  volume  loss 
(Skoog  et  al,  1985;  Koyama  et  al,  1988,  1992;  Malpas  et  al, 
1998).  With  progression  of  shock,  a  decrease  in  RSNA  may 
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occur  (Skoog  et  al,  1985;  Koyama  et  al,  1988, 1992)  as  a  sign 
of  irreversible  shock  (Koyama  et  al,  1988). 

Because  direct  tonic  nerve  activity  measurements  are 
invasive,  they  are  not  usable  as  clinical  monitoring  tools. 
Alternatively,  analysis  of  the  R-to-R-interval  (RRI)  by  fast- 
Fourier  transform  has  been  proposed  as  an  indirect  index  of 
cardiac  autonomic  modulation  (Akselrod  et  al.,  1981;  Pagani 
et  al.,  1986).  The  low-frequency  component  of  the  RRI 
power  spectrum,  or  low-frequency  (LF)  power,  is  influenced 
by  both  sympathetic  and  parasympathetic  activity.  High- 
frequency  power  (HF),  by  contrast,  is  influenced  by 
parasympathetic  activity  alone.  The  ratio  of  LF  to  HF  (LF/ 
HF)  has  been  proposed  as  an  index  of  “sympatho-vagal 
balance”  (Pagani  et  al.,  1986).  The  concept  of  sympatho- 
vagal  balance,  however,  has  been  critically  appraised  (Parati 
et  al.,  1995;  Eckberg,  1997;  Notarius  and  Floras,  2001).  At 
present,  the  indirect  heart  rate  variability  (HRV)  metrics  are 
said  to  reflect  cardiac  autonomic  influences  in  a  qualitative 
manner,  and  their  full  meaning  remains  controversial  (Parati 
et  al.,  2006). 

Despite  this  controversy,  HRV  metrics  are  widely  used  for 
obtaining  indirect  non-invasive  insight  into  cardiovascular 
autonomic  modulation  (Anonymous,  1996).  Association  of 
HRV  metrics  with  volume  status  (Butler  et  al.,  1994; 
Palazzolo  et  al.,  1998;  Goldstein  et  al.,  1999;  Batchinsky 
et  al.,  2007)  and  with  mortality  in  ICU  (Winchell  and  Hoyt, 
1996)  and  trauma  patients  (Cooke  et  al.,  2006)  has  suggested 
their  potential  utility  as  clinical  monitoring  tools  (Winchell 
and  Hoyt,  1996;  Grossman  et  al.,  2004;  Cooke  et  al.,  2006; 
Parati  et  al.,  2006). 

Saul  et  al.  (1990)  reported  a  positive  correlation  between 
the  LF/HF  ratio  and  directly  measured  sympathetic  nerve 
activity  in  humans,  while  others  have  not  (Kingwell  et  al., 
1994;  Pagani  et  al.,  1997;  Cooke  et  al.,  1999;  Cooke  and 
Dowlyn,  2000).  Kingwell  et  al.  (1994)  reported  dissociation 
between  the  LF  and  cardiac  sympathetic  activity  evaluated 
by  norepinephrine  spillover.  In  animal  studies  that  did  not 
involve  direct  RSNA  recordings,  RRI  LF  increased  during 
hemorrhage  in  conscious  dogs  (Madwed  and  Cohen,  1991), 
decreased  during  hypovolemic  shock  in  anesthetized  rabbits 
(Goldstein  et  al.,  1999),  and  showed  a  differential  response 
in  tachycardic  and  bradycardic  hemorrhaged  dogs  (Kawase 
et  al.,  2000). 

In  this  study  we  sought  to  characterize  the  changes  in 
RSNA  and  its  morphology  occurring  in  response  to  severe 
HS  in  sheep,  and  to  determine  the  relationship,  if  any, 
between  the  indirectly  assessed  LF/HF  ratio  of  the  RRI  and 
directly  measured  tonic  RSNA.  We  hypothesized  that  HS 
would  cause  concurrent  increases  in  both  the  LF/HF  ratio 
and  RSNA,  possibly  followed  by  withdrawal  of  RSNA  in 
pre-terminal  subjects.  We  also  hypothesized  that  HS  would 
be  associated  with  loss  of  RRI  complexity,  a  new  metric 
derived  from  non-linear  biosignal  analysis  methods,  which 
quantifies  RRI  irregularity  by  assessment  of  structural 
changes  within  the  signal  (Pincus,  1995;  Porta  et  al.,  1998; 
Richman  and  Moorman,  2000). 


2.  Materials  and  methods 

This  study  was  approved  by  the  U.S.  Army  Institute  of 
Surgical  Research  Animal  Care  and  Use  Committee  and  was 
carried  out  in  accordance  with  the  guidelines  set  forth  by  the 
Animal  Welfare  Act  and  other  federal  statutes  and  regulations 
relating  to  animals  and  studies  involving  animals  and  by  the 
1 996  Guide  for  the  Care  and  Use  of  Laboratory  Animals  of  the 
National  Research  Council. 

2.1.  Animal  preparation 

Neutered  male  sheep  (w  =  8),  weighing  35.9 ±  1.9  kg  were 
fasted  for  24  h  prior  to  the  study,  premedicated  by  ketamine, 
induced  with  isoflurane  via  mask  and  endotracheally 
intubated.  With  the  animal  in  dorsal  recumbency,  surgical 
anesthesia  was  carried  out  with  a  combination  of  isoflurane 
(0.47 ±0.03%)  and  intravenous  ketamine  (468 ±3.7  pg/kg/h) 
(KetaVed,  Ketamine  HCL,  Phoenix,  Scientific,  Inc.,  St. 
Joseph,  MO).  Animals  were  ventilated  in  volume-control 
mode  with  a  tidal  volume  of  13  ml/kg,  respiratory  rate  of 
12  per  min  and  fraction  of  inspired  oxygen  of  100%.  Tidal 
volume  and  respiratory  rate  were  adjusted  to  maintain  the 
end  tidal  level  of  CO2  between  30  and  40  mm  Hg.  Bilateral 
surgical  catheterizations  of  the  right  external  jugular  vein  and 
both  femoral  arteries  and  veins  were  performed  followed  by 
a  midline  laparotomy  and  splenectomy  to  exclude  the 
possibility  of  splenic  auto-resuscitation  during  hemorrhage. 
Each  animal  then  received  a  1-1  bolus  of  lactated  Ringer’s 
(LR)  solution  followed  by  a  maintenance  infusion  of  LR  at 
2  ml/kg/h.  Next,  the  left  kidney  was  exposed,  the  renal  nerve 
dissected  and  a  unipolar  Tungsten  metal  microelectrode  with 
a  75 -pm  shaft  (FHC,  Bowdoinham,  ME)  inserted  into  the 
nerve.  After  the  bursts  were  visually  verified  the  electrode 
was  glued  in  place  with  cyanoacrylate  adhesive  (Vetbond, 
3M,  Saint  Paul,  MN). 

2.2.  Physiologic  measurements 

Arterial  plasma  lactate  was  measured  with  an  i-STAT 
blood  gas  analyzer  (Abbot  Laboratories,  East  Windsor,  NJ). 
Standard  four-lead  EKG  was  obtained  from  surface  electro¬ 
des  glued  to  the  previously  depilated  skin  in  the  axillae  and 
outer  thighs.  A  continuous  cardiac  output  pulmonary  arterial 
catheter  (Model  746HF8,  Edwards  Lifesciences  LLC,  Irvine, 
CA)  was  placed  via  the  left  external  jugular  vein.  Arterial 
pressure  was  transduced  (Transpac®  IV,  Abbott,  North 
Chicago,  IL)  at  heart  level.  All  physiologic  waveforms  were 
displayed  continuously  using  a  clinical  monitor  (Viridia 
CMS  2000,  Boebingen,  Germany). 

2.3.  Experimental  protocol 

After  a  steady  nerve  signal  was  obtained  (verified  visually 
and  by  auscultation)  the  experimental  protocol  was  started. 
Animals  were  bled  in  four  10-min  cycles,  each  followed  by  a 
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15-min  observation  period.  Blood  was  manually  withdrawn 
via  a  syringe  from  the  femoral  arterial  line  during  each  1 0-min 
bleeding  period  at  a  constant  rate  of  1  ml/kg/min  for  a  total 
bleed  volume  (all  four  cycles)  of  40  ml/kg.  The  removed  blood 
was  collected  into  commercial  blood-collection  bags  contain¬ 
ing  CPD  anticoagulant  (Temmo  Corp.,  Tokyo,  Japan).  After 
the  fourth  hemorrhage  segment  and  observation  period,  sheep 
were  resuscitated  with  3  times  the  shed  blood  volume  with  LR 
over  a  period  of  25  min,  were  observed  for  15  min,  and  were 
given  their  shed  blood  back  via  infusion  over  25  min. 
Resuscitation  was  accomplished  with  the  Model  245  Ranger 
(Blood/Fluid)  Warming  system  (Augustine  Medical  Inc,  Eden 
Prairie,  MN).  Following  a  final  15-min  period  after  re- 
inftision,  surviving  animals  were  euthanized  by  means  of  an 
overdose  of  sodium  pentobarbital  (Fatal-Plus,  Dearborn,  MI). 

2.4.  Nerve  signal  analysis 

RSNA  signals  were  amplified  (90,000  times)  and  filtered 
at  low-  and  high-frequency  cut-offs  of  0.7  and  2  kHz 
respectively  using  a  Nerve  Traffic  Analysis  System  (Model 
662C-4;  University  of  Iowa  Bioengineering,  Iowa  City,  lA). 
The  presence  of  bursts  was  verified  both  visually  and  by 
auscultation.  Automatic  amplitude-based  detection  of  sym¬ 
pathetic  bursts  was  performed  with  WinCPRS  software 
(Absolute  Aliens  Oy,  Turku,  Finland).  Off-line,  the  signal  in 
all  files  was  brought  to  a  common  baseline  via  global 
polynomial  fitting  with  the  order  of  4  and  the  trends  were 
removed  using  the  Butterforth  high  pass  filter  with  the  cut¬ 
off  frequency  of  0.005.  RSNA  was  calculated  with  an 
averaging  window  width  of  5  s  and  average  window  step  of 
1  s  using  mean  normalization.  Discrete  RSNA  data  files  were 
selected  during  the  observation  periods,  and  were  normal¬ 
ized  with  the  baseline  activity  as  the  reference  value.  A 
scaling  factor  was  applied,  calculated  as  the  ratio  of  the  mean 
area  under  the  curve  for  one  burst  in  a  given  time  point  data 
file,  to  the  mean  area  under  one  burst  at  baseline  as  shown  in 
the  formula  below: 

SF  =  (area  TP/count  TP) /(area  BL/count  BE), 

where  SF  is  scaling  factor,  area  is  area  under  each  burst, 
count  is  burst  count  in  the  analyzed  segment,  TP  is  time  point 
of  experiment,  and  BE  is  baseline.  Next,  time  point  activity 
was  multiplied  by  the  scaling  factor.  RSNA  is  reported  as 
arbitrary  units. 

2.5.  Waveform  analysis 

EKG  and  arterial  blood  pressure  waveforms  were  con¬ 
tinuously  monitored  throughout  the  experiment  and  recorded 
at  500  Hz  using  an  in-house  digital  data  acquisition  system 
during  15-min-long  observation  periods  that  followed  each 
intervention  (i.e.  episodes  of  blood  withdrawal,  resuscitation 
or  blood  re-infiision).  This  allowed  us  to  do  the  analysis  on 
the  most  stable  discrete  data  sets.  Thus  7  time  points 
(baseline,  after  bleed  1 ,  after  bleed  2,  after  bleed  3,  after  bleed 


4,  after  ER  resuscitation  and  after  blood  re-infiision)  were 
analyzed  for  every  subject.  Next,  off-line,  800-beat  data 
sections  from  within  each  observation  period  were  imported 
into  the  WinCPRS  software.  Automatic  identification  of  R- 
waves  was  carried  out  by  means  of  an  iso-electric  line-shift 
algorithm.  We  verified  correct  identification  of  each  R-wave 
in  every  data  set.  The  analyzed  RRI  data  was  free  of  ectopy. 
The  software  generated  the  instantaneous  RRI  and  SAP  time 
series.  Before  spectral  calculations,  signals  were  windowed 
with  Hanning  function  and  re-sampled  with  a  5  samples/ 
second  rate.  The  re-sampled  signals  were  zero-padded  to  the 
next  power  of  2  before  the  FFT  routine  was  used  to  calculate 
the  power  spectra.  Next,  the  low-frequency  (LF)  power 
(frequency  range:  0.04-0.15  Hz),  and  high-frequency  (HF) 
power  (frequency  range:  0.15-0.4  Hz),  components  of  the 
RRI  were  calculated  for  each  data  set  with  baseline 
normalization.  Because  of  the  large  inter-animal  variability 
in  absolute  values  of  the  spectral  data  we  normalized  the 
spectral  data  to  the  total  power  (Anonymous,  1996).  All  data 
reported  in  this  study  are  normalized  spectral  data.  Normal¬ 
ization  was  achieved  by  dividing  the  integrated  LF  and  HF 
spectra  by  the  total  power  (minus  oscillations  occurring 
below  0.05  Hz)  and  multiplying  this  value  by  100.  The 
normalized,  integrated  areas  within  the  LF  band  were  divided 
by  normalized  areas  in  the  HF  band  to  derive  the  normalized 
ratio  (LFnu/HFnu).  This  ratio  was  taken  as  an  index  of 
sympatho-vagal  balance  (Pagani  et  al.,  1986). 

For  non-linear  analysis,  approximate  entropy  (ApEn), 
sample  entropy  (SampEn)  and  fractal  dimension  by  curve 
lengths  (FDCL)  were  calculated  for  the  same  800-beat 
segments  using  the  software.  Eight  hundred  beats  were 
selected  as  an  optimal  data  size  for  analysis  in  order  to  prevent 
any  bias  in  ApEn  which  is  prone  to  changes  in  the  number  of 
data  points  (Pincus,  1995).  Before  entropy  calculations, 
linear  trends  were  removed  in  all  segments  of  the  analyzed 
data.  For  both  ApEn  and  SampEn  calculations  the  dimension 
parameter  m  was  2  and  the  filter  parameter  r  was  20%  of  the 
standard  deviation.  Analysis  algorithms  are  identical  to  those 
used  by  Kuusela  et  al.  (2002). 

2.6.  Catecholamine  assays 

At  each  time  point,  5  ml  of  arterial  blood  was  collected 
into  vacutainer  tubes  containing  EDTA,  and  was  centrifuged 
at  1000  g  for  10  min.  The  obtained  plasma  was  stored  in 
polypropylene  tubes  at  -80  °C  until  the  day  of  the  assay. 
Plasma  samples  were  then  thawed  in  the  refrigerator  and 
assayed  according  to  He  et  al.  (1997)  with  modifications. 
Epinine  (25  pi  of  a  100  ng/ml  solution)  was  added  to  each 
1  ml  plasma  sample  and  served  as  the  internal  standard. 
Next,  1  ml  of  TRIS  buffer  (pH  9)  and  50  mg  of  acid-washed 
alumina  was  added,  followed  by  shaking  for  15  min.  Then, 
the  samples  were  centrifuged  at  1000  g  for  5  min  to 
precipitate  the  alumina,  which  was  washed  twice  with  1  ml 
of  TRIS  buffer.  The  alumina  was  mixed  with  150  pi  of  0. 1  M 
perchloric  acid  solution  and  vortexed  for  1  min,  than  the 
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Table  1 

Hemodynamic  and  catecholamine  data 


Time  point 
Event 

N 

0 

Baseline 

8 

1 

10  ml/kg 

8 

2 

20  ml/kg 

8 

3 

30  ml/kg 

8 

4 

40  ml/kg 

5 

5 

LR 

4 

6 

Re-infusion 

2 

(0  vs.  2) 

(0  vs.  4) 

(2  vs.  4) 

HR 

86±5 

82±3 

85±2 

117±12 

138±12 

114±11 

109±17 

NS 

<.0001 

<.0001 

SAP 

101±6 

74±5 

48±5 

38±3 

29±3 

99±18 

109±34 

<.0001 

<.0001 

<.0001 

DAP 

85±8 

59±6 

32±5 

19±3 

14±2 

59±13 

72±28 

<.0001 

<.0001 

<.0001 

CO 

2.7±0.1 

1.8±0.1 

1.2±0.1 

0.8±0.1 

0.5±0.1 

7.6±2.2 

6.8±0.7 

<.0001 

<.0001 

<.0001 

Lactate 

1.5±0.3 

1.6±0.3 

3.0±0.5 

7.5±1.1 

11.1±0.8 

15.6±1.1 

14.4±1.3 

.041 

<.001 

<.001 

Dop 

74±37 

85±47 

53±12 

118±23 

265±77 

164±31 

119±53 

NS 

.008 

.003 

Epi 

64±36 

45±8 

118±23 

3652±137 

14595±6149 

611±222 

262  ±106 

NS 

.002 

.002 

Norepi 

382±97 

429±88 

501±63 

1383±290 

6901  ±3052 

1095±210 

1298±374 

NS 

.009 

.011 

Time  point:  baseline  (0),  after  withdrawal  of  10  (1),  20  (2),  30  (3)  and  40  ml/kg  (4)  of  blood;  after  resuscitation  with  lactated  Ringer’s  (5);  after  re-infusion  of  shed 
blood  (6).  p^,  significance  level  when  baseline  compared  with  time  point  2;  p^,  significance  level  when  baseline  compared  with  time  point  4;  p^,  significance 
level  when  time  point  2  compared  with  time  point  4.  N,  sample  size  at  each  time  point;  HR,  heart  rate;  SAP,  systolic  arterial  pressure,  mm  Hg;  DAP,  diastolic 
arterial  pressure,  mm  Hg;  CO,  cardiac  output,  1/min;  lactate,  plasma  lactate  levels  in  mmol/1;  Dop,  plasma  dopamine,  ng/ml;  Epi,  plasma  epinephrine,  ng/ml; 
Norepi,  plasma  norepinephrine,  ng/ml.  Data  are  means  ±SEM. 


sample  was  centrifuged  at  1 000  g  for  5  min.  The  supernatant 
was  placed  in  microfilterfuge  tubes  and  centrifuged  for 
2  min,  and  100  pi  of  the  sample  was  injected  into  the  HPLC. 
For  calculation  of  norepinephrine,  epinephrine,  and  dopa¬ 
mine  levels  in  plasma  samples,  ratios  of  peak  areas  of  each 
catecholamine  to  the  peak  area  of  epinine  were  compared 
against  those  of  a  series  of  calibrators.  The  concentrations  of 
the  catecholamines  were  expressed  in  ng/ml.  The  HPLC 
system  consisted  of  a  Waters  515  pump.  Waters  717 
autosampler,  Coulochem  2  electrochemical  detector 
equipped  with  a  5011  detector  cell,  and  an  Alltima  Cl 8 
column  (5  pm,  4.6  (ID)  x  150  (length)  mm).  The  mobile 
phase  was  75  mM  phosphate,  1.7  mM  octanesulfonic  acid, 
25  pM  EDTA,  and  5%  (v/v)  acetonitrile  in  Milli-Q  water. 
The  flow  rate  of  the  mobile  phase  was  1.5  ml/min. 

2.7.  Statistics 

Data  analyses  were  carried  out  using  SAS,  version  8.1 
(SAS  Institute  Inc.  1999).  Variables  measured  at  different 
times  were  analyzed  using  a  mixed  model  ANOVA  allowing 
for  time  as  a  fixed  effect  and  replicate  subject  as  a  random 


effect.  The  covariance  structure  for  the  mixed  model  was 
determined  using  the  Bayesian  Information  Criterion. 
Hochberg’s  step-up  Bonferroni  method  was  used  for  p 
value  adjustment  for  multiple  comparisons  among  time 
points.  Log  transformation  was  used  when  appropriate. 
Statistical  significance  was  determined  as p<.05. 

3.  Results 

Hemorrhage  resulted  in  severe  shock,  such  that  only  25% 
of  animals  survived  until  the  end  of  the  experiment.  All 
animals  survived  removal  of  30  ml/kg  of  blood  at  time  point 
3.  Three  animals  died  after  bleed  3;  one  animal  died  after 
bleed  4;  and  two  animals  died  after  LR  resuscitation.  The 
physiologic  data  are  summarized  in  Table  1 .  Heart  rate  did 
not  change  between  baseline  (86  ±  5  bpm)  and  after  bleeds  1 
(82±30)  or  2  (85±2),  but  peaked  at  138±12  after  bleed  4 
(/?<.0001).  SAP  decreased  from  101  ±6  mm  Hg  to  48 ±5 
after  bleed  2  (/?<  .0001)  and  reached  its  nadir  after  bleed  4  at 
29 ±3  mm  Hg  (p<.0001).  Diastolic  arterial  pressure 
decreased  progressively.  Cardiac  output  decreased  from 
2.7 ±0.1  1/min  at  baseline,  to  0.5  1/min  after  bleed  4 


Table  2 


Renal  sympathetic  nerve  activity,  linear  and  non-linear  analysis  data 


Time  point 
Event 

N 

0 

Baseline 

8 

1 

10  ml/kg 

8 

2 

20  ml/kg 

8 

3 

30  ml/kg 

8 

4 

40  ml/kg 

5 

5 

LR 

4 

6 

Re-infusion 

2 

p' 

(0  vs.  2) 

(0  vs.  4) 

(2  vs.  4) 

RSNA 

586±61 

848  ±163 

1060±201 

1014±312 

216±91 

1866±787 

1435±1305 

NS 

.006 

.001 

LFnu 

0.59±0.08 

0.80±0.05 

0.88±0.04 

0.80±0.10 

0.87±0.06 

0.40±0.06 

0.57±0.13 

.009 

.001 

NS 

HFnu 

0.35±0.07 

0.17±0.03 

0.11±0.03 

0.16±0.07 

0.12±0.05 

0.52±0.06 

0.36±0.08 

.003 

.004 

NS 

LFnu/HFnu 

3.25±1.37 

6.57±1.74 

11.99±2.42 

9.31±1.72 

20.07±8.50 

0.83±0.17 

1.75±0.77 

.006 

.007 

NS 

ApEn 

0.91±0.11 

0.84±0.13 

0.82±0.14 

0.62±0.11 

0.59±0.09 

0.90±0.16 

0.99±0.03 

NS 

.033 

NS 

SampEn 

0.90±0.14 

0.85±0.17 

0.79±0.15 

0.57±0.11 

0.52±0.09 

0.83±0.22 

0.96±0.04 

NS 

NS 

NS 

FDCL 

1.72±0.04 

1.70±0.03 

1.65±0.03 

1.50±0.05 

1.49±0.06 

1.53±0.02 

1.75±0.06 

NS 

<.001 

.003 

Time  point:  baseline  (0),  after  withdrawal  of  1 0  (1 ),  20  (2),  30  (3)  and  40  ml/kg  (4)  of  blood;  after  resuscitation  with  lactated  Ringer’s  (5);  after  re-infusion  of  shed 
blood  (6).  p^,  significance  level  when  baseline  compared  with  time  point  2;  p^,  significance  level  when  baseline  compared  with  time  point  4;  p^,  significance 
level  when  time  point  2  compared  with  time  point  4.  N,  sample  size  at  each  time  point;  RSNA,  renal  sympathetic  nerve  activity;  LFnu,  normalized  low-frequency 
power  of  the  RRI;  HFnu,  normalized  high-frequency  power  of  the  RRI;  LFnu/HFnu,  ratio  of  the  normalized  low-frequency  power  to  normalized  high-frequency 
power  of  the  RRI;  ApEn,  approximate  entropy;  SampEn,  sample  entropy;  FDCL,  RRI  fractal  dimension  by  curve  lengths.  Data  are  means ±SEM. 
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Fig.  1.  Example  of  RSNA  activity  over  time:  A.  baseline;  B.  during  first  bleed  (500  ml  of  blood  removed).  Note  a  somewhat  increased  burst  amplitude  and  area 
when  compared  to  A;  C.  after  2nd  bleed  (980  ml  of  blood  removed).  Area  under  each  burst  increases  further.  Hyperpolarization  and  repetitive  firings  are  evident; 
D.  during  the  3rd  bleed  (1150  ml  of  blood  removed).  Formation  of  hyperpolarized  burst  conglomerates;  E.  after  the  4th  and  final  bleed  (1304  ml  of  blood 
removed).  RSNA  significantly  decreased;  F.  6  min  after  completion  of  4th  bleed.  Sympathetic  withdrawal;  G.  during  resuscitation  (1900  ml  LR  have  been 
infused).  Partial  return  of  RSNA;  H.  after  completion  of  resuscitation  with  LR,  3450  ml  infused).  Return  of  tonic  RSNA. 


(P  <.0001).  Lactate  increased  significantly  from  1.5 ±0.3  at 
baseline,  to  11.1  ±0.8  mmol/1  after  bleed  4  (/?<.001). 
Changes  in  base  excess  were  similar  (data  not  shown). 

Plasma  catecholamines  increased  significantly  after  bleed 
4  (Table  1). 

3.1.  RSNA  morphology 

RSNA  increased  in  7/8  animals  during  HS,  although  the 
time  point  at  which  the  highest  value  was  reached  varied 
among  animals.  Thus,  it  increased  numerically  but  not 
significantly  from  a  baseline  value  of  586±61  to  a  peak 
value  of  1060  ±201  after  withdrawal  of  20  ml/kg  of  blood 
(Table  2).  RSNA  then  decreased  in  all  animals,  reaching  a 
nadir  of  216  ±91  after  bleed  4  when  compared  to  both  post 
bleed  2  (/?  =  .006),  and  to  the  baseline  (/?  =  .001)  (Table  2). 


Fig.  2.  Changes  in  the  normalized  RRI  high-frequency  power.  *,  /7  =  .003, 
baseline  vs.  post  bleed  2.  f,  p  =  .0QA,  baseline  vs.  post  bleed  4.  See  text  for 
details. 


Fig.  1  shows  the  morphology  of  RSNA  at  each  time  point 
in  one  of  the  survivors  of  the  hill  protocol.  The  pattern  of 
multiple  burst  complexes  depicted  in  Fig.  3D  is  one  of  the 
most  pronounced  of  all  seen  in  the  study  and  was  present  in 
the  two  animals  that  survived  the  entire  study.  Similar 
changes — featuring  grouping  of  consecutive  bursts  into 
doublets  (Fig.  3B),  triplets,  and  multiple  bursts  (Fig.  3C) — 
were  observed  as  a  function  of  progressive  volume  depletion 
in  all  animals.  In  the  4  animals  surviving  past  the  LR 
resuscitation,  activity  increased  to  more  than  3  times  the 
baseline  value  (1866  ±787)  and  at  the  end  of  the  study,  in  the 
eventual  2  full  survivors  of  the  protocol  the  RSNA  remained 
almost  300%  of  the  baseline  value. 

3.2.  Spectral  analysis  of  EKG 

RRI  LFnu  increased  significantly  from  a  baseline  value  of 
0.59 ±0.08  to  0.88 ±0.04  after  bleed  2  (/?  =  .009)  and  was 


Fig.  3.  Changes  in  approximate  entropy.  1,7? <-033,  baseline  vs.  post  bleed 
4.  See  text  for  details. 
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sustained  at  a  higher  level  after  bleed  4  when  compared  to 
baseline  (Table  2).  RRI  HFnu  decreased  from  baseline  values 
(0.35±0.07)  to  a  nadir  (0.11  ±0.03)  after  bleed  2  (/?=.003, 
Table  2,  Fig.  2).  After  bleed  4,  HFnu  remained  low  compared 
to  baseline  (0.35 ±0.07  at  baseline  vs.  0.11  ±0.03  after  bleed 
2,  /?  =  .004)  but  was  not  different  from  the  nadir  point  after 
bleed  2  (0.12 ±0.05,  NS)  (Fig.  2).  The  LFnu/HFnu  ratio 
increased  significantly  after  bleed  2  (p  =  .006)  and  remained 
elevated  throughout  the  HS  period  and  significantly  higher 
than  baseline  after  bleed  4  (p=.007)  (Table  2).  With  the 
LFnu/HFnu  remaining  elevated,  and  RSNA  showing 
withdrawal  after  bleed  4,  a  significant  relationship  between 
the  LFnu/HFnu  ratio  and  RSNA  using  all  data  points  was  not 
revealed  by  linear  regression  (p  =  A2l,  =  0.048,  n  =  5l). 

3.3.  Non-linear  analysis  of  EKG 

Both  RRI  ApEn  and  SampEn  (Table  2)  followed  similar 
trends  in  this  study,  decreasing  numerically  after  each  blood 
withdrawal  and  reaching  the  lowest  values  after  bleed  4  at 
which  point  changes  in  ApEn  were  significant  (ApEn:  0.91  ± 
0.11  at  baseline  vs.  0.59±0.09,/?=.033,  Table  2,  Fig.  3).  The 
fractal  dimension  by  curve  lengths  (FDCL)  progressively 
decreased  with  each  blood  withdrawal  from  a  baseline  of  1 .72  ± 
0.04  to  a  nadir  of  1.49±0.06  following  bleed  4  (/?<  .001).  This 
FDCL  value  after  bleed  4  was  also  significantly  lower  when 
compared  to  post  bleed  2  (/?= .003)  (Table  2). 

4.  Discussion 

We  examined  changes  in  RSNA  in  response  to  a  highly 
lethal  model  of  hemorrhagic  shock  (HS)  in  anesthetized  sheep, 
and  assessed  the  relationship  between  non-invasive  indexes  of 
cardiac  autonomic  regulation  and  RSNA.  There  were  4 
principal  findings:  (1)  moderate  blood  loss  of  20  ml/kg  caused 
hypotension,  morphologic  evidence  of  RSNA  activation,  and 
increases  in  LFnu,  LFnu/HFnu  ratio  and  plasma  catecholamine 
levels;  (2)  severe  blood  loss  of  40  ml/kg  led  to  RSNA 
withdrawal,  sustained  LFnu  and  LFnu/HFnu  ratio,  and 
continued  increases  in  catecholamine  levels;  (3)  all  stages  of 
blood  withdrawal  were  associated  with  a  decrease  in  the  RRI 
HFnu;  and  (4)  blood  withdrawal  led  to  decreases  in  RRI 
complexity,  measured  by  ApEn  and  FDCL. 

RSNA  morphology  changed  as  a  fiinction  of  volume  status. 
Despite  some  degree  of  inter-animal  variability  in  RSNA  res¬ 
ponse  to  hemorrhage,  7  out  of  8  animals  demonstrated 
increasing  burst  frequency  with  volume  depletion.  The  ampli¬ 
tude  and  area  under  each  burst  also  increased  and  grouping  of 
bursts  into  doublets  and  triplets  took  place  as  more  blood  was 
withdrawn.  Eventually,  burst  conglomerates  formed  before  the 
onset  of  sympathetic  withdrawal  in  survivors.  To  our  know¬ 
ledge,  this  is  the  first  report  on  RSNA  burst  conglomerate 
formation  during  severe  hemorrhagic  shock.  Conglomeration  of 
bursts  may  indicate  the  accumulation  of  multiple  action  poten¬ 
tial  firings  (Tang  et  al.,  2003),  possibly  providing  for  a  stronger 
sympathetic  discharge  in  the  face  of  progressive  volume  loss. 


With  repetitive  and  superimposed  firings  hyperpolarization  may 
then  take  place,  which  may  lead  to  neuronal  hyperexcitability 
(Kandel  et  al.,  1991)  and  the  development  of  multiple  burst 
complexes  (Fig.  3C-D).  Development  of  these  complexes  may 
be  a  “last  resort”  compensatory  mechanism  during  pre-terminal 
shock,  and  may  be  caused  by  altered  neuronal  recruitment, 
altered  firing  patterns,  and/or  prolonged  neurotransmitter  release 
eventually  leading  to  receptor  saturation.  Such  multiple  burst 
complexes,  seen  in  the  two  full  survivors,  may  characterize 
subjects  with  adequate  autonomic  reserves  to  withstand  near- 
lethal  hypotension. 

Sympathetic  withdrawal  occurred  in  all  surviving  subjects 
after  blood  loss  of  40  ml/kg.  Skoog  et  al.  (1985)  found  short- 
lasting  sympathetic  activation  followed  by  powerfiil  sympa¬ 
thetic  inhibition  as  assessed  by  RSNA  measurements  in 
anesthetized  hemorrhaged  rats.  The  authors  proposed  that  the 
observed  sympathetic  inhibition  was  mediated  by  mechanically 
sensitive  cardiac  vagal  afferents  triggered  by  markedly 
decreased  cardiac  filling  (Skoog  et  al.,  1985).  Koyama  et  al. 
(1988)  reported  decreased  sympathetic  activity  in  their  studies 
on  anesthetized  dogs  subjected  to  hemorrhage,  and  proposed 
that  prolonged  brain  ischemia  causes  depression  of  sympathetic 
outflow  in  shock.  Sympathetic  withdrawal  was  suggested  to 
represent  a  fundamental  mechanism  for  the  development  of 
circulatory  shock  in  humans  subjected  to  lower-body  negative 
pressure  (Convertino  et  al.,  2004).  Our  findings  support 
observations  in  the  literature  that  severe  hemorrhage  may 
cause  a  failure  of  compensation  characterized  by  a  decrease  in 
sympathetic  nerve  activity. 

A  decrease  in  RRI  HFnu  occurred  during  HS,  consistent 
with  withdrawal  of  vagal  modulation  to  the  heart  (Palazzolo 
et  al.,  1998).  Previously,  we  identified  a  decrease  in  RRI  HF 
in  anesthetized  swine  subjected  to  isovolemic  hemorrhagic 
shock  (Batchinsky  et  al.,  2007).  Similarly,  Palazzolo  et  al. 
(1998)  documented  loss  of  vagally  mediated  RRI  oscilla¬ 
tions  during  hypotension  in  conscious  dogs. 

RSNA  and  LFnu/HFnu  ratio  did  not  correlate  in  our  study. 
This  dissociation  between  direct  microneurographic  measure¬ 
ments  of  sympathetic  tone  and  indirect  non-invasive  estimates 
of  autonomic  modulation  derived  from  heart  rate  variability 
analysis  has  been  reported  in  other  settings  as  well.  For  example, 
Pagani  et  al.  found  a  good  correlation  between  LFnu  of  HRV and 
LFnu  of  muscle  sympathetic  nerve  activity  (MSNA)  variability, 
but  not  between  LFnu  (or  LFnu/HFnu)  of  HRV  and  MSNA 
burst  rate  in  humans  (Pagani  et  al.,  1997).  Floras  et  al.  (2001)  did 
not  find  within-subject  correlations  between  MSNA  and  LF/HF 
ratio  during  lower-body  negative  pressure  (Floras  et  al.,  2001). 
Cooke  and  Dowlyn  (2000)  found  that  RRI  spectral  analysis  did 
not  adequately  reveal  subtle  changes  in  sympathetic  traffic 
during  head-down  tilt  as  assessed  by  muscle  sympathetic  nerve 
activity.  As  shown  by  Kingwell  et  al.  (1994)  in  patients  with 
heart  failure,  cardiac  sympathetic  nerve  firing  measured  with 
standard  spillover  techniques  and  LF  changed  in  opposite 
directions.  Thus,  the  dissociation  between  direct  tonic  and 
indirect  HRV  measures  of  autonomic  activity  is  a  predominant 
finding  in  the  literature.  Additionally,  the  concept  of  differential 
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activation  of  the  autonomic  nervous  system  in  various  vascular 
beds  may  be  an  additional  reason  for  the  discrepancy  between 
the  indirect  measures  of  autonomic  modulation  of  the  heart 
(HRV)  and  direct  measures  of  sympathetic  tone  (Ninomiya 
et  al,  1971;  Barman  et  al,  1984).  For  example,  Barman  et  al. 
(1984)  demonstrated  that  in  cats  the  reticular  formation  controls 
the  sympathetic  discharge  to  cardiac,  renal  and  carotid  nerves  in 
a  non-uniform  fashion.  Direct  cardiac  sympathetic  nerve  activity 
measurements  (Jardine  et  al.,  2005)  will  likely  provide  more 
decisive  evidence  on  the  relationship  between  indirect  HRV 
metrics  and  tonic  sympathetic  traffic  in  the  heart.  Meanwhile, 
however,  our  findings  support  the  evidence  in  the  literature  that 
changes  in  surrogate  measures  of  cardiac  autonomic  modulation 
such  as  LFnu/HFnu  ratio,  do  not  correlate  with  direct  measures 
of  tonic  sympathetic  discharge  in  the  renal  vascular  bed. 

In  recent  years,  non-linear  statistical  methods  have  in¬ 
creasingly  been  applied  to  the  analysis  of  the  EKG.  These 
methods  are  designed  to  analyze  dynamically  changing  non¬ 
stationary  data  (Anonymous,  1996)  such  as  those  explored  in 
this  study.  Complexity  is  a  fundamental  characteristic  of  normal 
physiologic  time  series,  such  as  the  RRI.  It  has  been  suggested 
that  loss  of  complexity  or  “decomplexification”  is  a  feature  of 
disease  (Goldberger  and  West,  1 987)  and  of  impaired  adaptation 
to  stress  (Lipsitz  and  Goldberger,  1992).  Complexity  can  be 
assessed  by  methods  such  as  approximate  entropy  (ApEn). 
ApEn  quantifies  the  probability  of  encountering  a  similar  pattern 
within  the  signal  at  the  next  incremental  comparison  (Pincus, 
1995);  it  is  affected  by  the  number  of  data  points.  SampEn  is  a 
similar  approach,  in  which  the  vector  comparison  with  itself  is 
removed;  this  makes  this  method  independent,  in  theory,  of  the 
number  of  data  points  (Richman  and  Moorman,  2000).  Entropy 
is  high  in  irregular,  complex  signals  and  decreases  as  the  signals 
become  more  regular.  In  our  study,  both  entropy  metrics  (ApEn 
and  SampEn)  showed  similar  trends  decreasing  numerically 
with  each  blood  withdrawal  and  reached  a  nadir  after  the  4th 
bleed  at  which  point  the  values  for  ApEn  were  significantly 
lower  when  compared  to  baseline  (changes  in  SampEn  were 
non-significant).  These  findings  are  consistent  with  our  previous 
work  that  identified  a  decrease  in  ApEn  and  SampEn  during 
isovolemic  hemorrhagic  shock  in  anesthetized  swine  (Batch¬ 
insky  et  al.,  2007).  Palazzolo  et  al.  (1998)  identified  a  decrease  in 
ApEn  during  hypotension  in  conscious  dogs.  Complexity  in  our 
current  study  changed  in  opposite  direction  to  the  indirect 
measures  of  autonomic  cardiac  modulation,  LFnu  and  LFnu/ 
HFnu.  This  is  in  agreement  with  Porta  et  al.  who  used  a  new 
entropy  calculation  approach,  the  corrected  conditional  entropy 
method,  and  found  loss  of  complexity  in  the  presence  of 
increased  LFnu  during  head-up  tilt  in  humans  (Porta  et  al.,  2000) 
and  during  nitropmsside  infiision  or  handgrip  (Porta  et  al., 
2007).  Corrected  conditional  entropy  is  a  promising  tool  which 
allows  for  analysis  of  shorter  data  sequences  without  a-priori 
selection  of  the  embedding  dimension  (Porta  et  al.,  1 998),  unlike 
ApEn  and  SampEn.  Further  studies  will  be  needed  to  define  the 
optimal  entropy  method. 

Another  measure  of  complexity,  FDCL,  is  a  fractal  image 
analysis  method  that  describes  changes  in  the  structure  of  the 


RRI  time  series  by  evaluating  the  number  of  short  segments  or 
“sticks”  necessary  to  follow  the  signal  curve  (Chau  et  al.,  1 993, 
Kuusela  et  al.,  2002).  A  higher  FDCL  value  indicates  an 
exponential  increase  in  the  number  of  sticks  required  to 
describe  the  structure  of  the  signal,  implying  higher  complex¬ 
ity.  The  decrease  in  fractal  dimension  metrics  identified  in  this 
study  has  been  previously  documented  by  us  during 
isovolemic  hemorrhagic  shock  in  swine  (Batchinsky  et  al., 
2007).  These  data  are  also  consistent  with  the  findings  by  West 
et  al.,  who  reported  a  decrease  in  the  fractal  dimension  of  the 
RRI  during  lower-body  negative  pressure  in  humans  (West 
et  al.,  2004). 

At  this  point,  loss  of  complexity  during  HS  maybe  related  to 
decreased  vagal  activity,  since  it  correlated  with  decreased  HFnu 
in  both  this  and  other  studies  (Palazzolo  et  al.,  1998;  Kuusela 
et  al.,  2002;  Batchinsky  et  al.,  2007).  Alternatively,  loss  of 
complexity  may  also  be  viewed  as  a  non-specific  marker  of 
physiologic  deterioration  under  conditions  of  hypovolemia  or 
critical  illness. 

5.  Conclusions 

In  an  anesthetized  ovine  model  of  highly  lethal  hemorrhagic 
shock,  moderate  hypotension  (blood  loss  of  20  ml/kg)  was 
associated  with  an  increase  in  the  RRI  LF/HF  ratio  (an  indirect 
frequency-domain  index  of  sympatho-vagal  balance),  and  a 
decrease  in  the  RRI  HF  (an  index  of  vagal  activity).  Directly 
measured  tonic  RSNA  featured  increases  in  the  burst  fre¬ 
quency,  amplitude,  and  area;  coupling  of  bursts  into  doublets 
and  triplets;  and,  finally,  formation  of  burst  conglomerates. 
Blood  loss  of  40  ml/kg  led  to  RSNA  withdrawal.  LF/HF  did  not 
correlate  linearly  with  RSNA.  RRI  complexity,  measured  by 
ApEn  and  fractal  dimension  methods,  decreased  after 
completion  of  the  40  ml/kg  hemorrhage.  Measures  of  RRI 
complexity  may  provide  useful  information  about  patient 
status  during  hypovolemia.  Further  work  will  be  needed, 
however,  to  understand  the  relationship  between  directly 
measured  autonomic  nerve  activity  and  various  indirect 
measures  of  autonomic  modulation  of  the  heart. 
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